Critical behavior of the planar magnet model in three dimensions 
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We study the critical behavior of the three-dimensional planar magnet 

model in which each spin is considered to have three components of which 
only the x and y components are coupled. We use a hybrid Monte Carlo al- 
gorithm in which a single-cluster update is combined with the over-relaxation 
and Metropolis spin re-orientation algorithm. Periodic boundary conditions 
were applied in all directions. We have calculated the fourth-order cumulant 
in finite size lattices using the single-histogram re-weighting method. Using 
finite-size scaling theory, we obtained the critical temperature which is very 
different from that of the usual XY model. At the critical temperature, we cal- 
culated the susceptibility and the magnetization on lattices of size up to 42 3 . 
Using finite-size scaling theory we accurately determine the critical exponents 
of the model and find that i/=0.670(7), 7/i/=1.9696(37), and (3/v=0. 515(2). 
Thus, we conclude that the model belongs to the same universality class with 

the XY model, as expected. 
64.60.Fr, 67.40.-w, 67.40.Kh 

I. INTRODUCTION 

Our understanding of critical phenomena has been significantly advanced with the devel- 
opment of the renormalization-group (RG) theory |l] . The RG theory predicts relationships 
between groups of exponents and that there is a universal behavior. In a second order 
phase transition, the correlation length £ diverges as the critical point is approached, and 
so the details of the microscopic Hamiltonian are unimportant for the critical behavior. All 
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members of a given universality class have identical critical behavior and critical exponents. 

The three-dimensional classical XY model is relevant to the critical behavior of many 
physical systems, such as superffuid 4 He, magnetic materials and the high-Tc superconduc- 
tors. In the pseudospin notation, this model is defined by the Hamiltonian 

H= -J 'EiSfSf + £%£%), (1) 

<ij> 

where the summation is over all nearest neighbor pairs of sites % and j on a simple cubic 
lattice. In this model one considers that the spin has two components, Si = (Sf,Sf) and 
Sf + Sf = 1. 

In this paper we wish to consider a three component local spin Si = (Sf, Sf, Sf) and the 
same Hamiltonian as given by Eq. (0) (namely, with no coupling between the z-components 
of the spins) in three dimensions. Even though the Hamiltonian is the same, namely, there 
is no coupling between the z-component of the spins, the constrain for each spin is (Sf) 2 + 
(Sf) 2 + (Sf) 2 = 1, which implies that the quantity (Sf) 2 + (Sf) 2 is fluctuating. In order to 
be distinguished from the usual XY model, the name planar magnet model will be adopted 
for this model. 

The reason for our desire to study this model is that it is related directly to the so-called 
model-F || used to study non-equilibrium phenomena in systems, such as superfluids, with 
a two-component order parameter and a conserved current. In the planar magnet model, 
the order parameter is not a constant of the motion. A constant of the motion is the z 
component of the magnetization. Thus, there is an important relationship between the 
order parameter and the z component of magnetization, which is expressed by a Poisson- 
bracket relation |J. This equation is crucial for the hydrodynamics and the critical dynamics 
of the system. One therefore needs to find out the critical properties of this model in order 
to study non-equilibrium properties of superfluids or other systems described by the model 
F. In future work, we shall use model-F to describe the dynamical critical phenomena of 
superfluid helium. Before such a project is undertaken, the static critical properties of the 
planar magnet model should be investigated accurately. 
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Although the static properties of the XY model with Si = (Sf, Sf) have been investi- 
gated by a variety of statistical-mechanical methods |3HTl~11, the system with Si = (S?, Sf, Sf) 



has been given much less attention. So far the critical behavior of this model has been stud- 
ied by high temperature expansion fT2j and Monte Carlo(MC) simulation methods Hl3,14 



High temperature expansion provides the value for the critical temperature and the critical 
exponents. In these recent MC calculations [ 13|, 14|| , only the critical temperature is deter- 



mined. These MC calculations were carried out on small size systems and thus only rough 
estimates are available. 

In this paper we study the three-dimensional planar magnet model using a hybrid Monte 
Carlo method (a combination of the cluster algorithm with over-relaxation and Metropolis 
spin re-orientation algorithm) in conjuction with single-histogram re-weighting technique 
and finite-size scaling. We calculate the fourth order cumulant, the magnetization, and the 
susceptibility (on cubic lattices L x L x L with L up to 42) and from their finite-size scaling 
behavior we determine the critical properties of the planar magnet model accurately. 

II. PHYSICAL QUANTITIES AND MONTE CARLO METHOD 

Let us first summarize the definitions of the observables that are calculated in our sim- 
ulation. The energy density of our model is given by 

< e >= E/V = \ £ < SfSf + SfS? >, (2) 

<ij> 

where V = L 3 and the angular brackets denote the thermal average. The fourth-order 
cumulant Ul{K) can be written as 

<r" Tn > 

where m = y(M^+My+M^) 1/2 is the magnetization per spin, M = Ei and K = J/(k B T) 
is the coupling, or the reduced inverse temperature in units of J. The fourth-order cumulant 
Ul(K) is one important quantity which we use to determine the critical coupling constant 
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K c . In the scaling region close to the critical coupling, the fourth-order cumulant Ul(K) as 
function of K for different values of L are lines which go through the same point. 
The magnetic susceptibility per spin \ is given by 

X = VK(< m 2 > - <m > 2 ), (4) 

where m is the magnetization vector per spin. 

The three-dimensional planar magnet model with ferromagnetic interactions J > has 
a second-order phase transition. In simulations of systems near a second-order phase transi- 
tion, a major difficulty arises which is known as critical slowing down. The critical slowing 
down can be reduced by using several techniques and what we found as optimal for our 
case was to use the hybrid Monte Carlo algorithm as described in Ref. [16|]. Equilibrium 



configurations were created using a hybrid Monte Carlo algorithm which combines cluster 
updates of in-plane spin components |T?| with Metropolis and over-relaxation of spin 



re-orientations. After each single-cluster update, two Metropolis and eight over-relaxation 
sweeps were performed [IE ]. The K dependence of the fourth-order cumulant Ul(K) was 



determined using the single-histogram re- weighting method JTP|. This method enables us 
to obtain accurate thermodynamic information over the entire scaling region using Monte 
Carlo simulations performed at only a few different values of K. We have performed Monte 
Carlo simulation on simple cubic lattices of size L x L x L with 6 < L < 42 using periodic 
boundary conditions applied in all directions and 10 6 MC steps. We carried out of the order 
of 10000 thermalization steps and of the order of 20000 measurements. After we estimated 
the critical coupling K c , we computed the magnetization and the magnetic susceptibility at 
the critical coupling K c . 

III. RESULTS AND DISCUSSION 

In this section, we first have to determine the critical coupling K c , and then to examine 
the static behavior around K c . Binder's fourth-order cumulant [T5|] Ul(K) is a convenient 
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quantity that we use in order to estimate the critical coupling K c and the correlation length 
exponent v. 

Near the critical coupling K c , the cumulant is expanded as 



T 



U L = U* + U 1 L 1/U (!-—) + ■ 



(5) 



Therefore, if we plot Ul(K) versus the coupling K for several different sizes L, it is expected 
that the curves for different values of L cross at the critical coupling K c . In order to find 
the K dependence of the fourth-order cumulant Ul(K), we performed simulations for each 
lattice size from L = 6 to L = 42 at K=0.6450 which is chosen to be close to previous 
estimates for the critical inverse temperature H0. The U L (K) curves were calculated 
from the histograms and are shown in Fig. [I] for L—12, 24, and 32. 
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FIG. 1. Fourth-order cumulant Ul{K) versus coupling K for lattice sizes L=12, 24, and 32. 

If one wishes to obtain higher accurary, then one needs to examine Fig. 1 more carefully 
and to see that the points where each pair of curves cross are slightly different for different 
pairs of lattices; in fact the points where the curves cross move slowly to lower couplings 
for larger system sizes. For the pair which corresponds to our largest lattice sizes L=32 and 
42, the point where they cross is K c m 0.64455. In order to extract more precise critical 
coupling K c from our data, we compare the curves of Ul for the two different lattice sizes 
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L and 11 = bL and then find the location of the intersection of two different curves Ul and 



Ujy. As a result of the residual corrections to the finite size scaling fIo| , the locations depend 
on the scale factor b = L'/L. We used the crossing points of the L—12, 14, and, 16 curves 
with all the other ones with higher V value respectively. Hence we need to extrapolate the 
results of this method for (lnfe) -1 — > using ({7&l/£/l)t=t c = 1- in Fig. we show the 
estimate for the critical temperature T c . Our final estimate for T c is 



T c = 1.5518(2), K c = 0.6444(1). 



(6) 



For comparison, the previous estimates are T c =1.54(l) [T3|jT4]] obtained using Monte Carlo 
simulation and T c =l. 552(3) [12| obtained using high-temperature series. The latter result 
obtained with an expansion is surprisingly close to ours. 
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FIG. 2. Estimates for T c plotted versus inverse logarithm of the scale factor h=L'/L. The 
extrapolation leads to an estimate of T c =\. 5518(2). 



In order to extract the critical exponent u, we performed finite-size scaling analysis of 
the slopes of Ul versus L near our estimated critical point K c . In the finite-size scaling 
region, the slope of the cumulant at K c varies with system size like L 1 ^, 

dU L 



dK 



L l/u . 



(7) 
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In Fig. |3] we show results of a finite-size scaling analysis for the slope of the cumulant. We 
obtained the value of the static exponent v, 



v = 0.670(7). 



(8) 



For comparison, the field theoretical estimate || is ^=0.669(2) and a recent experimental 
measurement gives u=0. 6705(6) |2"0|j . 
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FIG. 3. Log-log plot of the slopes of U near the crossing point versus L. The slope gives an 
estimate for the critical exponent ^=0.670(7). 

In order to obtain the value of the exponent ratio j/v, we calculated the magnetic 
susceptibility per spin x & t the critical coupling K c . The finite-size behavior for x a t the 
critical point is 



X 



L7 /V_ 



(9) 



Fig. ^ displays the finite-size scaling of the susceptibility x calculated at i^ c =0.6444. ^From 
the log- log plot we obtained the value of the exponent ratio j/u, 



7/1/ = 1.9696(37). 

^From the hyperscaling relation, dv = 7 + 2/3, we get the exponent ratio (3/v, 

P/ v = 0.515(2). 
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FIG. 4. Log-log plot of the susceptibility versus the lattice size L at the critical coupling 
if c =0.6444. The slope gives an estimate for the critical exponent 7 /V=l. 9696(37). 

The equilibrium magnetization m at K c should obey the relation 

m ~ L~ p ' v (12) 

for sufficiently larger L. In Fig. ^ we show the results of a finite-size scaling analysis for the 
magnetization m. We obtain the value of the exponent ratio (3/v, 

P/u = 0.515(2). (13) 

This result agrees very closely to that of Eq. (|i~l~l) obtained from the susceptibility and the 
fourth-order cumulant. 

In conclusion, we determined the critical temperature and the exponents of the planar 
magnet model with three-component spins using a high-precision MC method, the single- 
histogram method, and the finite-size scaling theory. Our simulation results for the critical 
coupling and for the critical exponents are K c =0. 6444(1), u=0. 670(7), 7/V=0. 9696(37), and 
j3/v=0. 515(2). Our calculated values for the critical temperature and critical exponents are 
significantly more accurate that those previously calculated. Comparison of our results with 
results of MC studies of the 3D XY model with two-component spins [[7|,p|-4Tl[] shows that 
both the system with Si = (Sf, Sf) and the planar magnet system with Si = (Sf, Sf, S?) 
belong to the same universality class. 
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FIG. 5. Log-log plot of the magnetization versus the lattice size L at the critical coupling 
if c =0.6444. The slope gives an estimate for the critical exponent (3/v=0. 515(3). 



L 


X 


m 


12 


82.39(28) 


0.26195(55) 


14 


111.88(36) 


0.24219(43) 


16 


145.12(59) 


0.22567(55) 


18 


182.91(52) 


0.21241(35) 


20 


224.08(85) 


0.20072(49) 


22 


272.23(60) 


0.19163(23) 


24 


322.35(98) 


0.18308(32) 


32 


571.0(4.0) 


0.15833(66) 


42 


972.0(4.8) 


0.13749(40) 



TABLE I. Results for the magnetization and the susceptibility 



9 



IV. ACKNOWLEDGEMENTS 

This work was supported by the National Aeronautics and Space Administration under 
grant no. NAG3-1841. 



[1] K.G. Wilson, Phys. Rev. B 4, 3174 (1971). 

[2] P.C. Hohenberg and B.I. Halperin, Rev. Mod. Phys. 49, 435 (1977). 
[3] J.C. Le Guillou and J. Zinn-Justin, Phys. Rev. B 21, 3976 (1980). 
[4] D.Z. Albert, Phys. Rev. B 25, 4810 (1982). 

[5] P. Butera, M. Comi, and A. J. Guttmann, Phys. Rev. B 48, 13987 (1993); R. G. Bowers and 
G. S. Joyce, Phys. Rev. Lett. 19, 630 (1967). 

[6] Y.-H.Li and S. Teitel, Phys. Rev. B 40, 9122 (1989). 

[7] A. P. Gottlob and M. Hasenbusch, Physica A 201, 593 (1993); A. P. Gottlob, M. Hasenbusch 
and S. Meyer, Nucl. Phys. B (Proc. Suppl.) 30, 838 (1993). 

[8] N. Schultka and E. Manousakis, Phys. Rev. B 51, 11712 (1995). 

[9] N. Schultka and E. Manousakis, Phys. Rev. B 52, 7528 (1995). 

[10] W. Janke, Phys. Lett. A 148, 306 (1990). 

[11] Martin Hasenbusch and Steffen Meyer, Phys. Lett. B 241, 238 (1990). 
[12] M. Ferer, M.A. Moore, and M. Wortis, Phys. Rev. B 8, 5205 (1973). 
[13] B.V. Costa, A.R. Pereira, and A.S.T. Pires, Phys. Rev. B 54, 3019 (1996). 
[14] S.K. Oh, C.N. Yoon, and J.S. Chung, Phys. Rev. B 56, 13677 (1997). 
[15] K. Binder, Z. Phys. B 43, 119 (1981). 

10 



[16] P. Peczak and D.P. Landau, Phys. Rev. B 47, 14260 (1963). 

[17] U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

[18] F.R. Brown and T.J. Woch, Phys. Rev. Lett. 58, 2394 (1987). 

[19] A.M. Ferrenberg and R.H. Swendsen, Phys. Rev. Lett. 61, 2635 (1988); ibid, 63, 1195 (1989). 

[20] L.S. Goldner, N. Mulders, and G. Ahlers, J. Low Temp. Phys. 93, 131 (1992). 



11 



